Signatures of Disoriented Chiral Condensates from Charged Pions 
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We show that the variance in the number of charged pions (in a suitable range of momentum 
space) provides a signature for the observation of a disoriented chiral condensate (DxC). The signal 
should be observable even if multiple domains of D%C form provided the average number of pions per 
domain is significantly large than unity. The variance of the number charged pions alone provides a 
signal which can be used even if the number of neutral pions cannot be measured in a given detector. 
If the neutrals can be measured, however, the fluctuations in the total number of pions provides 
a signature which distinguishes disoriented chiral condensates from other hypothetical sources of 
coherent states of pions. 



I. INTRODUCTION 



During the past several years, there has been considerable excitement about the possibility of the formation of 
disoriented chiral condensates (D%C) in heavy ion collisions [[[] ^0|. The basic scenario is as follows: In an ultrarela- 
tivistic heavy ion collision some region thermalizes at a temperature above the chiral restoration temperature. If the 
system cools sufficiently rapidly back through the transition temperature, the region will remain in a chiral restored 
phase. However, this phase is unstable; small fluctuations in any chiral direction (a, ff) will grow exponentially. This 
can create regions where the pion field has macroscopic occupation. It should be stressed that this scenario is not 
j derivable directly from the underlying theory of QCD and contains a number of untested dynamical assumptions, 
^{■^ principally that the cooling is rapid. Thus the failure of the system to form a D^C cannot be used to rule out that 
the system has reached the chiral restoration temperature. On the other hand, observation of the formation of a D%C 
Q\ , would be clear evidence that the phase transition had been reached. 

Unfortunately, since the scenario is not derived directly from a well-defined theory, it is difficult to know precisely 
what constitutes observation of a DxC. Assuming the system forms a single large domain of D%C containing a large 
■ number of pions there should be clear signatures. In the first place, one expects an excess in the number of low px 
i— H | pions produced. They would be at low px since, by hypothesis, the region is large so the characteristic momentum 
is small; the excess would be measured relative to a purely statistical thermal distribution. Such a signal has the 
advantage of working even if multiple regions of D\C form provided each region is large enough so the characteristic 
momentum is sufficiently small to provide a discernible signal over the thermal background. Such a signal is not 
decisive since one could imagine some other collective low energy effects which produce low px pions. 

A much stronger signal of a single large domain of D%C has been proposed. Since the pions formed in a D^C 
are essentially classical they form a coherent state. The coherent state has some orientation in isospin space (or 
more precisely the system is a quantum superposition of coherent states with different orientations and particular 
correlations to the isospin of the remainder of the system In essence all of the pions in the domain are 

pointing in the same isospin direction. Provided the total number of pions in the domain is large, this implies that 
the distribution of the ratio of neutral to total pions in the domain is given by @-§||ir||,|l8| 

m = m m 

where R is the ratio of the number of 7To's in the D%C divided by the total number of pions and P(R) is the 
probability. The derivation of P(R) is quite simple and will be discussed below. The distribution in Eq. ([!]) is 
qualitatively distinct from a purely statistical distribution in which the emission of charged and neutral pions is 
uncorrelated. The distribution from uncorrelated emissions in the infinite particle number limit approaches a delta 
function at R = 1/3. For finite (but large) particle number the statistical distribution is narrowly peaked about 1/3 
with a variance, (R 2 ) — (R) 2 = ^ where Af is the total number of pions. Since these two distributions are so radically 
different one should in principle have a very clear signal if a single region of D%C where to form in heavy ion reactions 
and if the pions from the D%C are kincmatically separated from other pions in the system. 

The dramatic nature of the preceding signature is based in large measure on the assumption that a single large 
domain of D%C is formed. A priori this seems rather unlikely for the following reason: If a large region of the 
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system starts in a hot chirally restored phase and then rapidly cools through the phase transition, then there will be 
a large region which is unstable against growth of the pion field. Presumably, this happens as a "seed" fluctuation 
in a small region which rapidly grows. It takes a time of at least L/c for information about the formation of the 
domain to propagate a distance L. However during the time this fluctuation is growing out to L, the pion field at L 
has been sitting in an unstable situation. The characteristic time it can remain in this unstable configuration is r, 
the exponential growth time. If the information about the initial seed does not reach L is a time comparable to r 
the region near L will likely begin its own exponential growth but in a chiral direction uncorrelated from the initial 
growth. Thus, one expects domains of characteristic size cr |l3| , ^4| . 

The effect of multiple domains on the R distribution is fairly clear: it will tend to wash out the signal. If a 
large number of domains form and the pions emerging from different domains cannot be distinguished kinematically 
it is clear from the central limit theorem that the R distribution will approach a normal distribution. This normal 
distribution may be distinguished from the normal distribution arising from uncorrelated emission; the case of multiple 
domains of D%C will have a substantially larger variance. 

Unfortunately, there is an important practical limitation which makes it difficult to exploit the R distribution as 
a signature. Even under the most optimistic of scenarios, the total number of pions coming from D^C's will be a 
small fraction of the total number of pions. If one includes all pions produced in the reaction, the signal from the 
pions from the D^C will presumably be overwhelmed. Thus, it is highly desirable to use kinematic consideration to 
enhance the contributions coming from the D^C. In particular, it is sensible to study the R distribution for a sample 
restricted to low px pions only. In any scenario where the D%C is well defined, i.e., the occupation number is large 
is likely to require a moderately large regions of D%C and the characteristic momentum spread in the D^C will be 
fixed by the inverse size of the region. Thus one expects D^Cs to preferentially produce moderately low px pions. 
(One also should restrict the pions in the distribution to a moderately narrow rapidity window). 

As an experimental matter, it should be relatively straightforward to cut on the momentum of the charged pions 
to select low px pions in a given rapidity window. For neutral pions, however, it is not a simple matter. The neutral 
pions will decay in flight and will ultimately be detected as photons. If one is simply interested in the overall R 
distribution, without cuts, and if the detected photons come predominately from 7r° decays then one can use the n 7 /2 
as a surrogate for n^o. Recent experimental searches have exploited this strategy pl|-|2^]. However, in order to study 
the R distribution in a limited kinematical region it is necessary to reconstruct the 7r^ momenta from the observed 
photons in order to make kinematical cuts on the 7r° momenta. Since the number of neutral pions per event is large, 
the reconstruction of neutral pion momenta is likely to be a formidable task. 

This raises the following interesting question: Can one find a signature for the presence of regions of D^C of 
essentially the same quality as the R distribution but which does not require the measurement of neutral pions? In 
this article, we will show that the distribution of the number of charged pions (in a kinematically limited region) 
contains essentially the same information about D^C formation as the the R distribution. This should greatly aid in 
searches for D%C formation. 

We also discuss additional information that can be inferred if 7r°'s can be reconstructed. In particular, we show 
that the distribution of the total number of pions (in a limited kinematic region) provides a means to distinguish D^C 
formation from other hypothetical mechanisms for the production of a pion coherent state. 



We begin by studying an overly simple model and in subsequent sections we will generalize our results to more 
realistic scenarios. In this simplified situation we assume that in every collision a single large domain of D%C is formed 
with a large particle number. Moreover, we will assume that the field strength and spatial distribution of this domain 
do not vary event by event, and that the pions produced in the D%C are kinematically completely distinguishable from 
all other pions in the system (including those pions produced from "ex's" — i.e., fluctuations in the (qq) directions). 
Finally, we will assume that both isospin violating effects and explicit chiral symmetry breaking are negligible. 

By hypothesis, the region of D^C contains many particles and is essentially classical in nature. To simplify discussion 
we will adopt the usual convention of describing the physics in terms of the degrees of freedom in a linear sigma model 
with 0(4) symmetry, i.e., a and ir rather than directly in terms of the QCD degrees of freedom. We wish to stress, 
however, that we are not relying on the detailed dynamics of any particular variant of the a model. 

Being a coherent state, the D^C can be written in the following form: 



II. A SIMPLIFIED MODEL 




(2) 
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where = (7r|, 7rJ, 7Tq, a^) is the vector of creation operators in the chiral space. The unit vector n — 

(sin-0 sin# cos <fi, smip sin sin0, sin ip cos 0, cos-0) denotes the orientation of the DxC in the chiral space, f(k) is the 
distribution of the D%C in momentum space, and |vac) is the vacuum. 
The number operators for neutral and charged pions are 

n Q = ttIttq, n± = irln x + (3) 

and one can easily find their expectation values , i.e., (^o,±)(V , j 0) = (DxC(t/>, 9, <^)|no,± |DxC(V>, 6, </>)). It turns out 
that («o,±) can be factorized into the following form. 



(no J± >foM,0) = (n>flb,±(lM,0), ( 4 ) 



where lo 2 = k 2 



The expectation value (n) measures the total number of 7r's produced by the D%C if fully oriented in a pionic direction. 
The geometrical factors go(ip, 9, 4>) an d 9±(4>, 9, (p) will be called the neutral and charged proportions respectively, and 
they take the following forms. 

go{ip, 9, 4>) — sin 2 ip cos 2 9, g± (ip, 9, cp) = sin 2 ip sin 2 9. (6) 

Note that (n) does not depend on the orientation angles (ip, 9, <ft) while go,± do not depend on the dynamical variables 
f(k) and u. 

One can also evaluate the higher moments of hq_±, which gives 

« ± >0A, 9,4>) = (n 2 ) gl ± W,e,<l>), (7) 

with 

(n 2 ) = (n) 2 + (n). (8) 
Let's define the deviance S[X] of a distribution of variable X such that 

(X 2 ) = (1 + 5[X})(X) 2 , or (AX) 2 = (X 2 ) - (X) 2 = S[X}(X) 2 . (9) 

In this case, 

S[n] = l/(n) -> when (n) -> oo. (10) 

Actually it is straightforward to show that n is described by a Poisson distribution, which always has a small deviance 
as the variance is proportional to the mean. 

The above analysis shows that, for each set of orientation angles (ip,9,(p), the distributions of no,± are normal. 
However, since the orientation is randomly generated in the process of spontaneous symmetry breaking, one cannot 
predict (ip, 9, </>). On the other hand, since we are neglecting explicit chiral symmetry breaking, the system is equally 
probable to point in any direction in chiral space. Moreover, since by hypothesis we are in a semiclassical situation 
(large (n)), it is legitimate to work with probabilities rather than amplitudes. Using the technology of Ref. Jlq], the 
probability distribution in the angular variables is given by the unit measure: 

d 3 P(ip,6,cp) = sin 2 V> sin0 dip d9dct>. (11) 
2tt z 

One can use Eq. (^) to reparametrize the probability distribution (|ll|) in terms of the neutral and charged proportions 
9o,± = ( n o,±) I ( n ) i *- e -i the fraction of neutral or charged pions among all particles produced by the D^C. 

d 2 P(g , g±) = - , = 1 = dg dg±. (12) 
77 V3o(l-.9o-.9±) 

From this one can obtain the marginal probability distribution of go and g± by integrating over the other variable. 
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dP(go) =fo(go)dg = -\ dg , 

Ti" V go 

dP(g±) = f±(g±)dg ± = dg±. (13) 

These distribution functions are plotted in Fig. 1 in solid curves. It is obvious the both distributions are far from being 
normal. The function /o is heavily skewed towards the low end and actually diverges as 1/^/go when go — > 0. On the 
other hand, f± is fiat, and g± is equally likely to assume any value between and 1. This is drastically different from 
pion emission from an uncorrelated source (the dotted curves in Fig. 1 R), where both distributions would be normal. 




FIG. 1. The probability distribution functions /o,±(<?o,±)- The plot on the left is /o(<7o), and the right is f±(g±)- The solid 
curves are for DCC emission, while the dotted curves are for independent emission with n = 50. 



More quantitatively, one can calculate the first and second moments of go and g± . 

(go) = 1/4, (gl) = 1/8, 

(g ± ) = 1/2, (gl) = 1/3, (14) 

and hence the respective deviance 5[go] and <5[g±], 

6\go] = 1, %±] = 1/3. (15) 

What does the distributions of the proportions <7o,± tell us about the distributions of rio,±? It is obvious that when 
ii is fixed, ga.± give the pion distribution. In reality, of course, n is not fixed; we have shown that it behaves like a 
Poisson distribution if we have a conventional Glauber coherent state |2^|. However, since Poisson distributions are 
sharply peaked, the dispersion of n will simply sm68J" the distribution of tiq.-k slightly, without changing the overall 
shape qualitatively. More specifically, note that |^| 

( n o,±) = (n) (go,±), 

«±> = (n 2 ) (gl,±) = (1 + *[n]) (1 + %o,±]) «n) (, 9 o,±)) 2 = (1 + %o,±])K ± } 2 ■ (16) 
The first equality gives 

(n ) = (n)/4, (n±) = (n)/2, (17) 



For uncorrelated emissions, the probability distributions are Poisson-Gaussian with mean 1/4 and 1/2 for neutral and 
charged pions respectively. The variances depend on the number of independently emitted pions; the plots correspond to the 
case of n = 50. 

2 The two set of brackets in the right-hand side of the equations below have different physical origins. The expectations 
of <?o,± are statistical in nature, while that of n is via quantum mechanical smearing of a coherent state. Consequently the 
distributions of n and go,± are assumed to be uncorrelated, and we derived the equalities below. 
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which have the simple interpretation that, given the symmetry between the four directions (7r x , n y , 7To, a) in the chiral 
space, a quarter of the particles produced by the D%C will be 7To, while half of them will be n x or 7r y . On the other 
hand, the second equality gives 

<%o,±] = S[n] + 6[go,±] + S[n]S[g ,±}- (18) 

With S[n] = l/(ri), 6[go] — 1 and <5[<?±], one has 

S[no] = l + 2/(n), *[n±] = 1/3 + 4/3<n). (19) 

As expected, when (n) is large, the deviances <5[ro 0i ±] approach <5[<7o,±]- Note that both <5[rio,±] are of order 1, in 
contrast to an uncorrelated emission which would have S = 1/ (n). In an ideal world, such enhancements of S's would 
indicate the existence of D^C. 

Lastly, it is also useful to reproduce the aforementioned probability distribution of R = no/(riQ + n±) = <?o/(5o + .9±) 
by reparametrizing distribution (|l2| ) in terms of R and g± — go + <7± (the subscript t stands for total) . 



d 2 P{g u R) = - J J* dg t dR, (20) 

7r v «(i - gt) 



with the marginal distributions 



dP{gt) = ft(gt)dg t = -./ 9t dg t , 
7T V 1 ~ 9t 

dP(R) = f R (R)d R = -^j=dR. (21) 

Note that while /o is drastically skewed towards the low end and f± is flat, ft is skewed towards the high end. This 
may sound counter-intuitive, but one must bear in mind that the emissions of neutral and charged pions are not 
independent and there is no contradiction. Also note that fa is exactly as predicted in Ref. (l] U|]-||,|l8| (cf. Eq. ([!])). 



III. MORE REALISTIC SCENARIOS 



We have studied the simple case of pion emission from a single huge (in the sense that (n) is large) D^C domain. 
As discussed in the introduction, this scenario is presumably not realistic. We will now proceed to study more realistic 
scenarios with multi-domain formation. The main point is, even though the probability distribution is smeared out 
because of the lack of alignment (in the chiral space) between the different domains, one feature survives, namely 
the large variance of the distributions. In particular, we will see that the variances for both neutral and charged 
productions are still much larger than that of an uncorrelated emission. 

Let's consider the case which we have N domains, each with the same (n). We are also assuming both (n) and N 
are much larger than unity. The total number of pions of each species produced is the sum of pions of that particular 
species produced in each domain, the distribution of which has been discussed in the previous section. 

N N 
i=l i=l 

By the central limit theorem, the probability distribution of Eno.± will approach normal distributions when N is 
large. However, we will see that the variances of the Gaussian distributions will be much larger for pion production 
from a D%C than those of uncorrelated pion emission. 

Since the pion production in each domain are independent, are independent random variables. Hence the 

mean of Sno,± is just the sum of the means of all ti!q ±, 

N 

(Eno,±) = X>£t> = N{n 0>± ) = Af ,±, (23) 

i=l 

(i) 

and the variance of the sum no,± is just the sum of the variances of each rig ± . 
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In other words, 



N 

(ASn , ± ) 2 - ^(An« ± ) 2 = 7V(An , ± ) 2 
i=l 

- iV(n , ± ) 2 <SK±] = (<y[no,±]/iV>A/J i± . (24) 
c rv i (ASn 0:± ) 2 % ,±] %o,±](n ,±) 

' [E ^ ±] = = — = aa . ± • (25) 

In comparison with uncorrelated pion production, with 5 = 1/JVo ± we see that for multi-domain D%C the deviances 
are enhanced by a factor of eo,± = <5[«o,±]("o,±)- 

6 = ((n) + 2)/4, c± = «n)+4)/6. (26) 

^4 priori (n) can take any value, but eo,± are larger than unity for any value of (n) > 2. Even for a very modest 
(n) = 8, eo = 2.5 and e± = 2, leading to substantial widening of the corresponding distributions, an observable 
signature of D%C formation. For larger values of (n), the broadening will be even more pronounced. 

While the above scenario describes D%C with multi-domains, a probable feature of D%C formation in the real world 
(if it happens at all), it is still unrealistic in assuming all the domains are of equal strength, i.e., with the same (n). 
Instead one expects (n) of different domains to fall under a certain probability distribution, which depends on the 
details of the model. Naturally, one questions if the signatures discussed above still survive under such circumstances. 

Let's consider the case with N domains, with different (nfi) ^> 1. Equation ( p3| ) becomes 

N 

(£n ,±) = £>o?±> = Njn^±) = JV ,±, (27) 



»=i 



where (no,±) is the average of (n l j_). Equation ( |24| ) becomes 

N N 



(ASn , ± ) 2 = ]T(Ar4%) 2 = Ev^L) 2 ^]- (28) 

i=l t=l 

By the inequalities <5[n ±] > 5[go,±] (cf. Eq.(|l9|)) and J2i ( n ) 2 < Ei^ 71 )) 2 /-^ (mean of squares is larger than square 
of mean), we have 



(A£n , ± ) 2 > (X>& 2 1 %0,±] > (l>&) ^.±1/^ = (^.9o,±]/^ 2 



' 2 ± . (29) 



In other words, 



, AEn ^) 2 %o,± % ,± (n ,±) , , 

o Sn ,± = -r= rj- > — — = T7 • (30) 

(Compare Eq. (|25|).) Again, the deviances are much larger than that of uncorrelated emission with S = l/7Vo,± when 
(n) ^> 1 by the following enhancement factors. 

eo > M/4, e± > M/6. (31) 

So we see that, even with domains of unequal strengths, the number of neutral or charged pions produced by a D%C 
will still have a much wider distribution than that from independent, uncorrelated emission. 

Lastly, one may also ask if the distribution of Eno,± will approach normal distributions when N is large in the 
case of domains with unequal strengths. In this case 's do not all fall under the same probability distribution and 
the most simple form of the central limit theorem does not apply. On the other hand, there are generalized forms 
of the central limit theorem, which state that as long as the probability distributions are sufficiently "well behaved" , 
the sum of iV random variables will still fall under a normal distribution when N — > oo. It is actually possible to 
argue that the distribution of En± does approach a normal distribution by the Lindeberg generalization of the central 
limit theorem. (See, for example, Sec. 6.E of Ref. p5[.) Whether the same conclusion holds for Sno is still an open 
question. 
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IV. DISCUSSION 



Let us recapitulate what we have shown in the previous sections. We have shown that one can calculate the 
probability distribution of the number of neutral or charged pions produced as a result of D%C formation. The 
resultant deviances 8's are not of the order 1/Af as in an uncorrelated emission, but are instead enhanced by factors 
e's which are of order (n). Seeing such enhancements of deviances would be signatures of coherent pion productions. 

One can understand the origin of such enhancements of statistical fluctuations of the number of neutral or charged 
pions by considering the following analogy. Consider two groups of gamblers playing roulette in a casino: N lawyers 
at $100 tables, and 1007V physicists at the $1 tables, where the odds are the same. If each lawyer and physicist is 
given n chips, of $100 and $1, respectively (so that the total amount given to the lawyers, JV = N x lOOn, is the 
same as that given to the physicists, Af = 1007V x n), and is required to bet all of them, the average loss will be the 
same for both groups as long as they are following the same betting strategies. However, it is easy to see that the 
statistical fluctuation of the loss of the lawyers would be much larger for that of the physicists. In other words, the 
amount of loss, as well as its standard deviation, is "quantized" in units of the value of the bets. The larger the bet, 
the larger the fluctuation. On the other hand, one can also turn the argument around; a discerning external observer 
can deduce, with the knowledge of the gambling strategies, the size of the bets of the lawyers from the statistical 
fluctuations of the lawyers' losses, and do likewise for the physicists as well. 

Just as the chips in a single bet of a lawyer share the same fate (either win or lose), all the pions in a single D%C 
domain share the same orientation in the chiral space. As a result, the fluctuation of the number of pions in each 
direction in the chiral space is enhanced by a factor proportonal to (n), the number of pions in each D%C domain. 
And by reverse argument, one can deduce whether coherent pion emission is taking place by measuring the fluctuation 
of the number of emitted pions. 

Our analysis has assumed that all the pions originate from coherent emissions, and each of these coherent states 
has (n^) 3> 1. In the real world, there is contamination from independent pion emissions, and the total numbers of 
neutral or charged pions are the sums of these two contributions. 

no,± — n Q ± + n a ± , (62) 

where "c" and "inc" stands for "coherent" and "incoherent" , respectively. The variances of the numbers of coherently 
produced neutral or charged pions are enhanced while those of incoherent production are not. 

(An^) 2 = e , ± <r4^>, (An«) 2 = (ng), (33) 

where eo,± are of order (n) ^S> 1 (cf. Eq. (|26|), (|3l|)). Then one can calculate the variance of the sum of the two 
contributions. 

(An) 2 = (Ar^i) 2 + (A^ n ± c) ) 2 = e , ± (n), e ,± = X* + (1 - x) = 1 + x(e - 1), (34) 

with 

X = (n£i>/<no,±>, 1 - X = (n^)/(n Q , ± ). (35) 

The parameter x measures the fraction of pions which are coherently produced: x = 1 when all the pions are from 
D^C, while x = when all of them are independently emitted. Obviously, the more incoherent pions in the sample, 
the smaller is the enhancement factor eo,±. 

While these incoherently emitted pions dilute our signatures for D%C, they have different momentum spectra from 
those from D^C. D%C pions, being produced from coherent state, carry low px- The typical pr is of the order of 1/i, 
where L is the size of the domain from which the pion originates. In contrast incoherently emitted pions can carry 
high pt- Therefore applying a low px cut can minimize the noise from incoherent pion emissions. 

It is also advantageous to measure the rapidity of the pions and count their numbers in narrow rapidity windows. 
Bear in mind that the rapidities of the pions are, up to small dispersions, equal to that of the original domain. As we 
have mentioned, it is probable that many domains are formed in a single collision, and all these domains may have 
different rapidities. For example, the domains at the surface of the "fire ball" are moving with high speed relative to 
the domains at the center. By binning the pions according to their rapidities, one can partially separate the pions 
from different domains, and the signals are enhanced as a result. 

In summary, we suggest the following procedure in looking for signatures of D^Cs. 

• Count the number of neutral or charged pions event by event from heavy ion collision experiments and measure 
their individual transverse momenta and rapidities. 



7 



• Apply a low px cut to suppress the noise due to uncorrelated pion emission. 

• Bin the events in different rapidity windows. 

• In each rapidity window, plot the number of events vs. the number of neutral or charged pions in histograms. 

• Evaluate the mean, (no,±), and the variance, (Ano^) 2 , in each rapidity window. 

• If we find (Ano^) 2 is substantially larger than (no,±), then we are seeing signatures from D^Cs. 

The above procedure allows us to search for signatures from D%Cs by counting only the charged pions. This is 
important as, with our present technology, it is difficult to count the number of 7ro's in a momentum bin, which would 
mean reconstructing all the pions from photons — a formidable task. On the other hand, with great experimental 
effort, it may be possible to count the neutral pions as well in the future. In that case, we will be able to distinguish 
D^C formation from other mechanisms of coherent pion productions. For example, one can count nt, the number 
of pions (both neutral and charged) in each rapidity window. For D^C formation, or any other mechanisms of 
coherent pion productions where the field is aligned with a random direction in the four-dimensional chiral space 
(7T x ,TT yi TTQ,a), the fluctuation of n t is large. On the other hand, for mechanisms of coherent pion productions where 
the field is aligned with a random direction in the three-dimensional isospace (jr x , ir y ,iro) but without involving the a 
direction, it is straightforward to show that the fluctuation of nt is small. 

In summary, we have constructed signatures for D^C formation in heavy ion collisions which do not require counting 
the number of neutral pions. Instead, we suggest counting the number of charged pions produced, and a large 
fluctuation would be a signal of D^C formation. We believe these new signatures will be useful in searches for D^C 
at RHIC and LHC. 

Support of this research by the U.S. Department of Energy under grant DE-FG02-93ER-40762 is gratefully ac- 
knowledged. 
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